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Abstract. We develop direct and inverse scattering theory for Jacobi op- 
erators with steplike coefficients which are asymptotically close to different 
finite-gap quasi-periodic coefficients on different sides. We give a complete 
characterization of the scattering data, which allow unique solvability of the 
inverse scattering problem in the class of perturbations with finite first mo- 
ment. 



1. Introduction 

In this paper we consider direct and inverse scattering theory for Jacobi oper- 
ators with steplike quasi-periodic finite-gap background, using the Marchenko [TS] 
approach. 

Scattering theory for Jacobi operators is a classical topic with a long tradition. 
Originally developed on an informal level by Case in [5], the first rigorous results 
for the case of a constant background were given by Guseinov [H] with further 
extensions by Teschl [TH] , [10] . The case of periodic backgrounds was completely 
solved in [24] (who in fact handle almost periodic operators with a homogenous 
Cantor type spectrum) respectively [8] using different approaches. Moreover, the 
case of a steplike situation, where the coefficients are asymptotically close to two 
different quasi-periodic finite-gap operators, was solved in [11] (see also pQ, [7]) 
under the restriction that the two background operators are isospectral. It is the 
purpose of the present paper to remove this restriction. 

We should also mention that scattering theory for Jacobi operators is directly 
applicable to the investigation of the Toda lattice with initial data in the above 
mentioned classes. See for example [3J, [Ej, [13] for steplike constant backgrounds, 
and [S], [IH]) [13], [H]) and [H] for periodic backgrounds. For further possible 
applications and additional references we refer to the discussion in [11] . 

Finally, let us give a brief overview of the remaining sections. After recalling 
some necessary facts on algebro-geometric quasi-periodic finite-gap operators in 
Section [2 we construct the transformation operators and investigate the properties 
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of the scattering data in Section [3l In Section |4] we derive the GePfand-Levitan- 
Marchenko equation and show that it uniquely determines the operator. In addi- 
tion, we formulate necessary conditions for the scattering data to uniquely deter- 
mine our Jacobi operator. Our final Section [5] shows that our necessary conditions 
for the scattering data are also sufficient. 

2. Step-like finite-band backgrounds 

First we need to recall some facts on quasi-periodic finite-band Jacobi operators 
which contain all periodic operators as a special case. We refer to 20, Chapter 9] 
and [8] for details. 

Let be two quasi-periodic finite-band Jacobi operators^ 

(2.1) Hff{n) = af(n)f(n + 1) + af(n - l)/(n - 1) + b±(n)f(n), f e £ 2 (Z), 
associated with the Riemann surface of the square root 

2ff±+l 

(2.2) P±{z) = - [] \/ z - E ?> E± <Ef <-.-<E 2g±+1 , 

3=0 

where g± £ N and is the standard root with branch cut along (— oo, 0). In fact, 

are uniquely determined by fixing a Dirichlet divisor X)f=i(/ i j t ; a f)> wnere 
Hj € [E 2 j-i,E 2 j] and er* e { — 1, 1}- The spectra of consist of g± + 1 bands 

g± 

(2.3) v±:=v(H±) = \J [E% , E% +1 ]. 

3=0 

We will identify the set C \ a(H^) with the upper sheet of the Riemann surface. 
The upper and lower sides of the cuts over the spectrum are denoted by cr u and tr 1 
and the symmetric points on these cuts by A u and A , that is, 

/(A lI )=lim/(A + ie), /(A z ) = lim/(A — ie), \E* ± . 

ej.0 eiO 

We will develop the scattering theory for the operator 

(2.4) Hf(n) = a(n-l)f(n-l)+b(n)f(n) + a(n)f(n + l), neZ, 

whose coefficients are asymptotically close to the coefficients of on the corre- 
sponding half-axes: 

(2.5) M(VW - + \Hn) - &±(n)|) < oo. 

n=0 

The special case H~ = has been exhaustively studied in [8j (see also [24]) and 
the case where H~ and are in the same isospectral class cr_ = er + was treated 
in Several results are straightforward generalizations, in such situations we 

will simply refer to [S], [TT] and only point out possible differences. 
Let ip^(z, n) be the Floquet solutions of the spectral equations 

(2.6) H±t/](ri) = zijj(n), zeC, 



Everywhere in this paper the sub or super index "+" (resp. "— ") refers to the background on 
the right (resp. left) half-axis. 
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that decay for z G C \ <7± as n — ► ±00. They are uniquely defined by the condition 
ip^(z,0) = 1, ip^(z,-) 6 £ 2 (Z±). The solution ip+(z,n) (resp. ip~(z,n)) coincides 
with the upper (resp. lower) branch of the Baker- Akhiezer functions of (resp. 
H~), see [2D]. The second solutions i(£(z, n) are given by the other branch of 

the Baker- Akhiezer functions and satisfy (z, •) € £ 2 (Z T ) as z 6 C \ o±. Their 
Wronskian is equal to 

(2-7) W^(i>f(z),^(z))=±^j-r, 

q q q P±{z) 

where 

(2.8) p ± (z)- ^ 

satisfy by our choice of the branch for the square root 

(2.9) Im(p±(A u )) > 0, Jm(p±(X 1 )) < 0, A e <r ± . 
In (12. 7p the following notation is used 

(2.10) W±JJ,g) := a±(n) (f(n)g(n+l) - f(n + l)g(n)) . 

Note that ip^(z,n), ■t/j^ : (z,n) have continuous limits as z — > A"' 1 S clj.' 1 \ <9o±, 
where 

dcr± = ...,£; 2g±+1 }, 
and they satisfy the symmetry property 



(2.H) Vf(A 1 ,»)=Vg(A»,n)=^(A n > n), A G <r±. 

The points [f£j,Oj) } 1 < J < form the divisors of poles of the Baker- 

*g±- 



Akhiezer functions. Correspondingly, the sets of Dirichlet eigenvalues {/if, ju^ } 



can be divided in three disjoint subsets 

M ± = {pf | pf € R\<7± is a pole of i/;±(z, 1)}, 
(2.12) M± = {fif | /x± e K\cr ± is a pole of ^(z,l)}, 

M ± = { M ±| M ±e5a ± }. 
In order to remove the singularities of ip^(z,n), ip^(z,n) we introduce 
<5 ± (z):= J] 



(2-13) 



1 3 • 



s ± (z):= n i*-4) n v*-^, 
s±(z)~ n ( z -^) n 

where J} = 1 if there are no multipliers, and set 

(2.14) 4>±(z,n) = 5±(z)ijf(z > n), $£(z,n) = $±(z)il)f(z,n). 

Lemma 2.1. The Floquet solutions ip^, xf)^ have the following properties: 
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(i) The functions ip^ (z , n) (resp. ^(z,n)) are holomorphic as functions of z 
in the domain C \ (a± U M± ) ( resp. C \ (o~± U M± ) ), take real values on the 
set R \ a±, and have simple poles at the points of the set M± (resp. M±). 
They are continuous up to the boundary cr^_ U a± except at the points in 
M± and satisfy the symmetry property (|2 . 1 1 1) . For E G M±, they satisfy 



1 



Moreover, the estimate 



%(z,n) = O 



1 



, z^EeM ± . 



(2.15) ^{z,n)-^{E,n)=0{yf7=E), E G da±, 
is valid. 

(ii) The following asymptotic expansions hold as z — > ±oo 

n— 1 i n— 1 -i 

(2.16) ^(«,n)=^(n* a*(j)) (l±-£* ^ + I) + 0(^)), 



j'=o 



where 



it /c?) 



/ n— 1 

II /(?'): n > n 0, 
3=no 

1, n = n , £* /(i) = 



11-1 



no — 1 



J=™o 



/ 71— 1 

j=no 

0, n = no, 

no — 1 

- E /0')> 71 < "0- 



11 /(i) \ n< no, 

V j=n 

(iii) The functions ip^ (A, n) form a complete orthogonal system on the spectrum 
with respect to the weight 



(2.17) 



(2.18) 



(2.19) 



duj ± (X) = —p±{X)dX, 

27T1 



namely 



where 



tp q (X,m)ip (X,n)dw±(X) = S(n,m), 



f(X)dX := / f(X u )dX - / f(X l )dX. 

J (TV J (j\ 

Here 5(n,m) = 1 if n = m and S(n,m) = else is the Kronecker delta. 
3. Scattering data 



Now let H be a steplike operator with coefficients a(n), b(n) satisfying (|2.5 
The two solutions ijj±(z, n) of the spectral equation 



(3.1) 



Hip = zip, z6C, 



which are asymptotically close to the Floquet solutions tpf(z, n) of the background 



equations (|2.6p as n — > ±oo, are called Jost solutions. They can be represented as 
(see 0) 

±00 

(3.2) ip±(z,n) = ^2 K ±{ n ^ m )i ) q i z > TO )> 
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where the functions K±(n, .) are real valued and satisfy the estimate 
(3.3) 



±00 



\K ± (n,m)\<C ± (n) £ (|a(j) - af(j)\ + W) ~ ±™ > ±n > 0. 

The functions C±(n) > decrease monotonically as n — > ±00. Moreover, we have 
K+(n+ l,n + l) 



(3.4) 



a(n) = a,q(n) 
a(n) = a~(n) 



""" X+(n,n) 
AT- (n, ra) 



K-{n + l,n + ly 

= + <(n)^±^ - <(n - 1) = ^ 

M * K+[n,n) H K+(n — l,n — 1) 

, _ _ K-(n, n — 1) if_(n + l,n) 

6(n) - b q (n) + a q (n- 1) T ) }- ^ - a q (n) 



K-(n,n) q v 'K_{n + l,n+ 1) ' 
which implies (cf. [8]) the following asymptotic behavior of the Jost solutions as 



z -> ±00 using 03), (GDI]), 
(3.5) 

n-1 



^ ± (z,n) = z^X ± (n,n)(n* (l+(s±(n)± &±(j-o)) +0 ( ? ; 

3=0 3=1 

where 



±00 



(3.6) B±(n)= (^M-6M). 



m— nil 

J 



For A € cr± U a second pair of solutions of (|3.1p is given by 

±00 

(3.7) #t(A,n) = ^±(n,m)^(A,m), Ae^U^, 



which cannot be continued to the complex plane. Note that ijj±(\,n) = ip±{A,n), 
A € a±, and from (|2.5|) . (|3.2[) we conclude 



(3.8) ^(^ ± (A),^±(A)) = iy±(^ ± (A),^ ± (A)) =±p±(A)- 1 . 

The Jost solutions ip± are holomorphic in the domains C\ (cr± U M±) and inherit 
almost all properties of their background counterparts listed in Lemma 12.11 As 
before, we set 

(3.9) 4>±(z,n) = 5±(z)ip±(z,n), rp±(z,n) = 5±(z)ip±(z,n). 
The following Lemma is proven in [8]. 

Lemma 3.1. The Jost solutions have the following properties. 

(i) For all n, the functions ip±(z, n) are holomorphic in the domain C \ (cr± U 
M± ) with respect to z and continuous up to the boundary (crj U o± ) \ do~± , 



where 

. \ U \ , \ I , \ , Hi 1 \ \ ■ 1 



(3.10) i>±(X u ,n) = ip±(A\n), A e (ct± U a l ± ) \ da 
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The functions ip±(z, n) are real valued for z 6 R\cr± and have simple poles 
at ix j E M± . Moreover, ip± are continuous up to the boundary a±U a±. 
(ii) At the band edges we have for A € a± l 



V>±(A,n)- V> ± (A,n) = o(l), Eeda±\M ± , 
( - > 1>±(\n)+il>±{\,n) = o( 7 £%), E € M±. 

Next, we introduce the sets 

(3.12) cr (2) := <T + ncr_, a { ± ] = clos ((7± \ f7 (2) ), (7 := (7+ U(7_, 

where tr is the (absolutely) continuous spectrum of H and o+ U er_ resp. cr^ 2 ) are 
the parts which are of multiplicity one resp. two. We will denote the interior of the 
spectrum by int(cr), that is, int(er) := o~ \ da. 

In addition to the continuous part, H has a finite number of eigenvalues situated 
in the gaps, ad — {Ai,...,A p } c M \ a (see, e.g., [18]). For every eigenvalue we 
introduce the corresponding norming constants 

(3.13) 7 ^ 1 fc = ^|^ ± (A fc ,7 1 )| 2 , l<k<p. 

Moreover, ip±(Xk,n) = c^ T (A fe ,n) with c^c^ = 1. 
Let 

(3.14) W{z) := W(i>-(z),ip+(z)) 

be the Wronskian of two Jost solutions. This function is meromorphic in the domain 
C \ a with possible poles at the points M + U M_ U (M+ n M_ ) and with possible 
square root singularities at the points M + U M_ \ (M + n M_). Set 

(3.15) W(z) = W$-(z),$+(z)), W{z) = W$-(z),j>+(z)), 

then W(A) is holomorphic in the domain C \ R and continuous up to the boundary. 
But unlike to W(z) and W(z), the function W(X) may not take real values on the 
set M \ a and complex conjugated values on the different sides of the spectrum. 
That is why it is more convenient to characterize the spectral properties of the 
operator H by means of the function W , which can have singularities at the points 
of the sets M + U Af_ . We will study the precise character of these singularities in 
Lemma 13.21 below. 

Note that outside the spectrum the function W{z) vanishes precisely at the 
eigenvalues. However, it might also vanish inside the spectrum at points in <9er_ U 
da + . We will call such points virtual levels of the operator H, 

(3.16) a v := {E e a : W(E) = 0}, 

and we will show that a v C da U (da^ n da_ ) in Lemma 13721 All other points E 
of the set da+ U <9er_ correspond to the generic case W(E) ^ 0. 

Our next aim is to derive the properties of the scattering matrix. Introduce the 
scattering relations 



ii.i 



(3.17) T T (X)ip ± (\,n)=^ T (\,n)+R T (X)i; T (\,x), A e , 

where the transmission and reflection coefficients arc defined as usual, 



(3.i8) t ±( a)._ w{MXU±wy R±W -~ wiMX )^ ±{ x)y 



SCATTERING THEORY FOR STEPLIKE QUASI-PERIODIC BACKGROUND 



The equalities in (|3.18|) imply the identity 
1 1 



W(X), AG(J (2 \ 



T + (X) P+ (X) T_(A)p_(A) 

where W(X) is the Wronskian of two Jost solutions (|3.14|1 . This Wronskian plays 
an important role in the characterization of the properties of the scattering matrix. 
Namely, the following result is valid. 

Lemma 3.2. The entries of the scattering matrix have the following properties: 
I. 



(a) T ± (A u )= r ± (A'), Aea±, 
i?±(A u )=i?±(A 1 ), Xea ± , 

(b) =# = * ± W, Ae4 1} , 
1 ±[X) 

(c) i-|i? ± (A)P = ^|r ± (A)| 2 , Aea( 2 \ 

(d) R±(X)T ± (X) + R T {X)T ± (X) = 0, X e . 

II. The functions T± (A) can be extended analytically to the domain C \ (crU M± U 
M±) and satisfy 

(3-19) 1 = 1 — = W(z). 

T+{z)p+{z) T-(z)p-(z) 

The function W(z) has the following properties: 

(a) The function W(z) = 6 + (z)5-(z)W(z) is holomorphic on C \ er with simple 
zeros at the eigenvalues Xk, where 

fdW ,\ 2 1 



(3-20) hH Afe )/ 

\ dz J 7+,fe7-,fe 

Moreover, 



(3.21) W(X U ) = W(X ), Aea, W(z) e E, zeR\a. 

(b) The function W (z) = S + (z)S^(z)W(z) is continuous on the set<C\o~ up to the 
boundary a u Ucr l . It can have zeros on the set daU(da^ Dda_ ) and does not vanish 
at the other points of the spectrum a. IfW(E) — as E £ da U (9crv n 9<r^), 
i/ien 

(3.22) -J— = O | . ) , where X e a close to E. 

(c) In addition, 

(3.23) T+(oo) =T_(oo) > 0. 

III. (a) The reflection coefficients R±(X) are continuous functions on int(cr^_' 1 ) . 
(b) If E € <9cr+ D <9cr_ and W^E 1 ) ^ 0, iften £/ie functions R±(X) are also 
continuous at E. Moreover, 



(3.24) R ± (E) 



-1 /or £ £ M±, 
1 for E e M±. 



s 
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Proof. I. The symmetry property (a) follows from formulas (I3.18|) and (13.10|) . For 



(b), use p.ISp and observe that ip=(\) are real valued for A G mt(<r±^). Let 



A € int(cr^). By (l3~T71) . 

\T±\ 2 W(i> T ^) = (\R ± \ 2 - l)W(iP ± M, 

and property (c) follows from (|3.8[) . The consistency condition (d) can be derived 
directly from definition (|3.18[) . 

II. The identity (|3T9|) follows from (|3T8| . (a) The Wronskian inherits the 
properties of ip±( z )i so it remains to show (|3.20|) . If W(z ) — for z E C \ er, then 

(3.25) 4>±(zo,n) = ^^(zo,^ 

for some constants (depending on zq), which satisfy c~c + = 1. In particular, 
each zero of W (or W) outside the continuous spectrum is a point of the discrete 
spectrum of H and vice versa. 

Let 7± j be the norming constants defined in (|3.13[) for some point of the discrete 
spectrum Xj. By virtue of [2U], Lemma 2.4, 

±W$.(z),i, + (z))\ x . =W n (^-(A i ),^+(A J -)) + W„(^-(A i ),^+(A i )) 

(3.26) =-Vi(A i) fc)^ + (A i( fc) = -^, 

Since cjc+ = 1, we obtain (|3.20|) . 

(b) Continuity of W up to the boundary follows from the corresponding property 
of ip±(z,n). We begin with the investigation of the possible zeros of this function 
on the spectrum. 

First let A € int(cr( 2 )) := <j( 2 ) \ da (2 \ that is, 5-^0 a nd 6+ j- 0. Sup- 
pose Fy(Ao) = 0, then ?/; + (Ao,n) = c-0_(Ao,n) and Tp+(Xo 1 n) = c-0_(Ao,n), i.e. 
W / (V'+,'0+) = \c\ 2 W(ip~,ijj-). But this implies opposite signs for p+,p~ by (|3.8p . 
signp + (A ) = — sign p_(A ), which contradicts (|2.9p . 

Let Ao G int((7j_ 1 ' ) ) and I^(Ao) = 0. The point Ao can coincide with a pole 
/i G Mzf. But ^±(^0;^) an d tp±{^0i n ) are linearly independent and bounded, 
and ?/> T (A ,n) G E. If W(X ) = 0, then tp T — cf = ip± which implies 
W(tp±,ip±)(Xo) — 0, a contradiction. 

In the general mutual location of the background spectra the case Xq = E E 
{do-W n int(cr±)) C int(cr) can occur. If W(E) = 0, then W(tp±,ip^){E) = 0, 
where ip T are defined by (|3.9[) . The values of ipzp(E, •) are either purely real or 
purely imaginary, therefore W(ip±, ip^)(E) = 0, that is, ip±(E,n) and tp±(E,n) are 
linearly dependent, which is impossible at inner points of the set o~±. 

Thus, the virtual level a v of H defined in 1|3.16[) can only be located on the set 
da- n da + , that is, 

(3.27) a v C da U (<9cr a) n da^) . 

To prove (|3.22p . take E E a v and assume for example E E <r+. By (|3.17p and 

(CT , 

S+(X)^(X,n) 



5-(X)p+(X)W(X) 



S+(X)i; + (X,n)+R + (X)^ + (X,n). 
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Choose no such that ip-(E,no) ^ 0. By continuity we also have ^>_(A, no) 7^ in 
a small vicinity of E. Then 



5+(X)tp + (X,n ) + R + (X)i p+(X,n ) = 



tp- (A, n 
Correspondingly, 



w(X) \5 2 + (x)^x~eJ \Vx^eJ' 

which proves |3~.22I 

(c) Equation (|3~23"|) follows from ((375]) . 

III. (a) follows from the corresponding properties of ip±(z) and from II, (b). To 
show III, (b) we use that by (|3.18p the reflection coefficients have the representation 



(3 - 28) * ±(A) = -w^im^m ~ T WW 

and are continuous on both sides of the set int(cr±) \ (Mip U M^). Moreover, 



|i?±(A)| 



W(X) 



where the denominator does not vanish on the set a± \ a v . Hence i?±(A) are 
continuous on this set since both numerator and denominator are. 

Next, let E e d<r± \ a v (in particular W(E) ^ 0). Then, if E £ M ± , we use 
(|3~2"gj) in the form 



(3.29) * ± (A) = -1 T M^M)). 

W(X) 



which shows R±{X) — ► —1 since ip±(X) — ip±(X) — > by Lemma 13.11 (2). This 
settles the first case in (]3 . 24[) . Similarly, if E <E M± , we use (|3.28l) in the form 



(3.30) R±[x) ^ l± ^m^W^ ±{ X)^ { X))^ 



which shows R±(X) — > —1 since S±(X) = 0(VA — E) and ip±(X) + ip±(X) = 
o(— =) by Lemma \3.1l (2). This settles the second case in (|3.24p as well. □ 

4. The Gel'fand-Levitan-Marchenko equation 

The aim of this section is to derive the inverse scattering problem equation (the 
Gel'fand-Levitan-Marchenko equation) and to discuss some additional properties 
of the scattering data which are consequences of this equation. 

Theorem 4.1. The inverse scattering problem (the GLM) equation has the form 

(4.1) K ± (n, m) + V K ± (n, l)F ± (l,m) = * ( "' m \ , ±m>±n, 
tL K ± {n,n) 
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where 

F±(m,n)=<b R±(X)ipf(X,m)ipf(X,n)dw± 



(4.2) + / \T T (\)\ 2 ^(\,m)^(\,n)dcJ T +Y,l±,k^(Xk,n)4>±(\ k ,m). 



=F fe=l 



Proof. Consider a closed contour T e consisting of a large circular arc and some 
contours inside this arc, which envelope the spectrum a at a small distance e from 
the spectrum. Let ±m > ±n. The residue theorem, (|2.17[) . (|3.5p . (|3.20[) . and 
equality $ T {X k ,n) = cJ^±(A fe ,n) imply 



1 / ^Mj^(M ,,_ <?(n,m) \p R / V> T (A,n)V>*(A,m) \ 
27ri/ re W(A) ~Jr±(n,n) ^ W(A) J 

(4.3) = P^^ -Y,7±,kMXk,n)i>f(\ k ,m), 

K±{n,n) y 

since the integrand is meromorphic on C\er with simple poles at the eigenvalues 
Afc and at oo if m = n. It is continuous till the boundary except at the points 
E e (9(7+ U da- where 

-0t(A, n)ip^(X,m) ( 1 \ 



by ()3.22j) . On the other hand, as e — * 0, 

1 /" <MA,n)-0^(A,m) 

— d\ — 



2iri J a W[X) 

1 I (V± (A, n) + i*± (A)^± (A, n)) ^ (A, m) 



27ri/ CT± T±(A)W(A) 
1 /" 4>^(X, n)il}„(\, m) 



-0± (A, n)^q (A, m)du ± + ® i?± (A)i/>± (A, ra)V>* (A, m)(iw ± 

J CT-J- 



dX. 



W{\) W(X) 
It remains to treat the last integrand. By (|3.17p and Lemma [3~^1 I, 



?V (A), 



^ (A, ») = T T (\)ip± (A, n) - R T (A)tA t (A, n) = T T (X)ip± (A, n) - =±=^ T (X, n), 

1-T\X) 

and therefore 



iyr T + iyr T t t t t 

— ! = — = — Wit ( n — ==w±(n) = — ==w± (n , 

W |VK| 2 T T TV ' W y ' W V 7 
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since WT T + WT T = 2Re(WT T ) = on a T . In summary, flO) and (03]) yield 

S(n,m) _ ^ , / D , xv .,. M „ X .,± M x, ,± 



K± (n, n) J a± q 

+ [ \T T (X)\^ ± (X,n)^(X,m)d^+j2l±J±(^^)^(X 1 ,m) 1 

"H 3=1 

and applying (|3.2p finishes the proof. □ 

As it is shown in [5], the estimate (|3.3p for K±(n,m) implies the following esti- 
mates for F±(n, m). 

Lemma 4.2. The kernel of the GLM equation satisfies the following properties. 
IV. There exist functions C±(n) > and q±(n) > 0, n £ Z±, such that C±(n) 
decay as n — > ±oo, |n|g(n) £ £ 1 (Z±), and 

±00 



\F±(n,m)\<C ± (n) ]T q(j), 

j=n-\-m 

±oo 

(4.6) H\F±(n,n) - F±(n ± 1, n ± 1)| < oo, 

n— no 

|n| |a^(n)_F± (n, n + 1) — a^(n — l)F±(n — 1, n)| < oo. 



In summary, we have obtained the following necessary conditions for the scat- 
tering data: 

Theorem 4.3. The scattering data 

S = \r+(X), T+(X), X £ cr^ 1 ; R_(X), T_(A), A £ a^ 1 ; 

(4.7) 

Ai,...,A p £ R\ (cr+ Ucr_), 7±,i,---,7±,p 

satisfy the properties I-III listed in Lemma \3.S\ The functions F±(n,m), defined 
in \4 .3\), satisfy property IV in Lemma\4~J^ 



In fact, the conditions on the scattering data given in Theorem 14.31 are both 
necessary and sufficient for the solution of the scattering problem in the class (|2.5p . 
The sufficiency of these conditions as well as the algorithm for the solution of the 
inverse problem will be discussed in the next section. 



5. The inverse scattering problem 

Let be two arbitrary quasi-periodic Jacobi operators associated with se- 
quences a^(n),6^(n) as introduced in Section [5] Let S be given scattering data 
with corresponding kernels F±(n,m) satisfying the necessary conditions of Theo- 
rem 14.31 

First we show that the GLM equations (|4.1[) can be solved for K ± (n, m) if 
F±(n,m) are given. 
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Lemma 5.1. Under condition IV, the GLM equations have unique real-valued 
solutions K±(n, •) 6 £ (n, ±00) satisfying the estimates 

±00 

(5.1) \K±(n,m)\ < C±(n) £ q(j), ±m > ±n. 



Here the functions q± (n) and C± (n) are of the same type as in l{4.6\ 
Moreover, the following estimates are valid 

±00 



(5.2) 

/ ±00 



\n\\K±(n,ri)-K±(n±l,n±l)\ < 00, 

i=n 
±00 

|n| \ a^(n)K± (n, n + 1) — a^(n — l)K±(n — 1, n)\ < 00. 



Proof. The solvability of (|4.ip under condition (|4.6[) and the estimates (|5.ip . (I5.2|) 
follow completely analogous to the corresponding result in [8] Theorem 7.5]. To 
prove uniqueness, first note that the GLM equations are generated by compact 
operators. Thus, it is sufficient to prove that the equation 



±00 



(5.3) f(m)+52F±(l,m)f(£)=0 

has only the trivial solution in the space ^(njioo). The proof is similar for the 
"+" and "— " cases, hence we give it only for the "+" case. Let f{£), £ > n, be a 
nontrivial solution of (|5.3p and set f(£) = for £ < n. Since F + (£, n) is real- valued, 
we can assume that f(£) is real- valued. Abbreviate by 

(5.4) /(A) = 53^+(A,m)/(m) 

the generalized Fourier transform, generated by the spectral decomposition (|2.18[) 
(cf. (2^). Recall that /(A) € L\ oc {a\ U er^)- 

Multiplying (|5.3p by /(m), summing over m £ Z, and applying (|2.18p , (|4.2p . 
(|5.4p . and condition I, (a), we have 



(5.5) 



U(A)|^ + (A) + 2Re / R + (X)f(\ydoj + (\) 
+ [ /(A) 2 |^(A)| 2 d^(A) + f] 7+ , fc f^^ + (A fc ,7 1 )/(n)) =0. 

•^ tT - k=l / 



The last two summands in ()5.5j) are nonnegative since /(A) 6 R for A 6 cr_ and 
"4>q{Xk) G IK- We estimate the first two integrands by 

|/(A)| 2 + Rei? + (A)/(A) 2 > |/(A)| 2 - |i? + (A)/(A) 2 | > (l - \R + {\)\) |/(A)| 2 
and drop the last summand in (|5.5p . thus obtaining 

(5.6) 2 / (1 - |i? + (A)|)|/(A)| 2 du, + (A) + / /(A) 2 |H(A)| 2 dcMA) < 0. 
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Here we also used that 

(l-|i?+(A)|)|/(A)| 2 ^ + (A)=0, 



(l),u 

+ 



which follows from condition I, (b). Since \R+(X)\ < 1 for A £ int(cr( 2 )) and 
u)-(X) > for A 6 int(er_ ) we conclude that 

/(A) = for A G cr (2) U cr a) = <7_ . 

The function f(z) can be defined by formula (|5.4[) as a meromorphic function on 
C \ <7+. By our analysis it is even meromorphic on C \ cr+ and vanishes on cr_. 
Thus f(z) and hence also f(m) are equal to zero. □ 

Next, define the sequences a±,b± by 

K + (n+ l,n+l) 



a+ (n) = a+(n)- 



a_(n) = a„ (n) 



K+(n,n) 
K— (n, n) 



(5.7) 



q y J K-(n+ l,n+ 1)' 



M n ) = + ngfo) ' \ a, (n - 1)- 



K+(n,n) q v #+(n- 1)' 
t / \ n . _/ ,X-(n,n-l) _ *f_(n + l,ri) 

* y iv_(n, n) y a_ (?i + 1, n + 1) 

and note that estimate (|5.2j) implies 

(5.8) n{\a ± (n) ~ af(n)\ + \b ± - b±(n)\} E t\Z±). 

Lemma 5.2. The functions ip±(z 1 n), defined by 

±00 

(5.9) ip±(z,n) = K±(n,m)i()g(z,m), 

m—n 

solve the equations 

(5.10) a±(n — l)ip±(z, n — 1) + b±(n)ip±(z, n) + a±(n)ip± (z,n+ 1) = zip±(z, n), 
where a±(n),b±(n) are defined by |5. 71 ) ■ 

Proof. Consider the two operator!! 

(H±y)(n) =a±(n - l)y±(n - 1) + b±(n)y±(n) + a±(n)y±(n + 1), n e Z. 
Define two discrete integral operators 

±oo 

(£±/)(n) = ^ tf ± (n,m)/(m). 

m—n 

Then (cf. [8]) the following identity is valid 

H±K.± = K± H q , 

which implies (|5. 10[) . □ 

The remaining problem is to show that a+(n) = a_(n), b + {n) = &_(n) under 
conditions II and III on the scattering data S. 



2 We don't know that H± is limit point at =poo yet, but this will not be used. 
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Theorem 5.3. Let the scattering data S, defined as in j4-ty j satisfy conditions I, 
(a)— (c), II, III, (a), and IV. Then each of the GLM equations j-l-lty has unique 
solutions K±{n,m), satisfying the estimate i5.2\) . The functions a±(n),b±(n) , de- 
fined by \5. 7| ), satisfy \5.8\) . 

Under the additional conditions III, (b) and I, (d), these functions coincide, 
a + {n) = a_(n) = a(n), b + {n) = fo_(n) = b(n), and the data S are the scattering 
data for the Jacobi operator associated with the sequences a(n),b(n). 



The proof of Theorem 15 . 31 takes up the remaining section and is split into several 
lemmas for the convenience of the reader. 

To prove uniqueness of the reconstructed potential we follow the method pro- 
posed in [TS]. Recall that, by Lemma 12.11 (iii), the functions ipg(X, n) form an 
orthonormal basis with corresponding generalized Fourier transform. Split the ker- 
nel of the GLM equation (|4.2[) into three summands F±(m,n) — F ry ±(m,n) + 
Fh.± (™, n) + Fd,±(m, n) and set 

(5.11) G±(n,m) :=J2 K ±( n > l )Fr,±(l,n). 

l—n 

Then one obtains as in [8l Theorem 8.2] that the functions /i=p(A, n), denned by 

( 5 - 12 ) T J mu \T T (£)\ ip±{^n) — du T (0 

p w± f.7,± f\.\ J,±( 



±2^ 7± > fc ^ ± ( Afc > n ) -frr — 

k=l fe 



satisfy 



(5.13) T ± (X)h T (X,n)=i/j±(X,n) + R±(X)ip±(X,n), X G 4 . 

Despite the fact that h T (X,n) are defined via the background solutions cor- 
responding to the opposite half-axis Z± , they share a series of properties with 
ip^f(X,n). Namely, we prove 

Lemma 5.4. Let h^{z,n) be defined by formula 15.12)) on the set o^' 1 . 

(i) The functions hz^(z,n) — S^(z)h^(z,n) admit analytic extensions to the 
domain C\o~. 

(ii) The functions hzf(z, n) are continuous up to the boundary tr"' 1 except pos- 
sibly at the points da + U <9<j_ . Furthermore, 

h zr(A u ,n) = h T (X\n) G R, XeR\a^, 

(5.14) . , 

h T (X n ,n) = h T (X l ,n), X G mt(cr T ). 

(iii) For large z the functions hz^(z,n) have the following asymptotic behavior 

(5.15) K±MT±{oQ) (jl a±(j)) (l + 0(-)), ,^oo. 
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(iv) We have 

W ± (h T (z), ip±(z)) ■=a±(n)(h T (z, n)tp±(z, n+ 1) — h T (z, n+ l)tp±(z, n)) 
= ±W(z), 

where W(z) is defined by i3.19\) . 

Remark 5.5. Note that we did not establish the connection between the function 
W(z) and the functions W ± (i/j+(z,n),4'^(z,n)), which can depend on n, because 
and ip- are the solutions of Jacobi equations corresponding to possibly different 
operators H + and H- . 

Proof, (i). To show that h T (z, n) have analytic extensions to C \ c, we study each 
term in (15.12|) separately. 

First of all, note that due to the representation 



1 5 T (z)^/UfJo +1 ^~E^ 

(5.16) T±(z) = 



p±{z)W(z) 

o±(z) W(z) 
the functions ^(z,n) = d T (z)( T (z,n) 1 where 

i)„{z, n) 

(5-17) £=(*,n):= 

can be continued analytically to C \ a. This also holds for the second term since 
G±(n,-) € € 1 (Z) are real-valued. 

Next we discuss the properties of the Cauchy-type integral in the representation 
(I5.12[) . We represent the third summand in (|5.12[) multiplied by T^ t 1 (z) as 

(5.18) e T (z,n):= T ^T / e T {z,^n) ^ 



where 



e T (z,Z,n) = (ffifc •). C=f(*, •)) 

(5-19) = IMOI 2 IM^ ^ ( ^ )w ± ( ^± (e ) c ( )}< 

p t (0\w(0\ 2 M£) 2 ' 

By property II, (a) the function W(£) has no zeros in the interior of <7^' u . Thus, 

for z , the functions 0^z(z, .,n) are bounded in the interior of cr^ and the 

only possible singularities can arise at the boundary. We claim 

iO(JT^E) for E£o- V: m 
(5.20) T (z,£,n) = l_Y \ ' , * Eeda¥,z^E. 



This follows from ^§f- = 0{^f£~~E) together with W{£) = 0(1) if E <£ a v and 



1/W(£) = (^(l/v 7 !^ 7 ^) by II, (b) if E G £7„. Therefore, T are integrable and the 
third summand of (|5.12[) also inherits the properties of C=f( z i n )- 

Finally, the last summand in (|5.12[) again inherits the properties of ^(z,n) 
except for possible additional poles at the eigenvalues A^. However, these cancel 
with the zeros of W(z) at z = 
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ii). We consider the boundary values next. The only nontrivial term is of course 



the Cauchy-type integral (|5.18[) as z — > A G 'mt(a^). First of all observe that by 
(12~TD and 1^15]) . 

where the functions S±W are bounded and nonzero for A G intfog 5 ) by II, (a). 
Hence the Plemelj formula applied to (|5. gives 

6 T (A, n) — ± T t — S : — «£, AGmt 4' ), 



where both terms are finite. Here j- denotes the principle value integral. Therefore, 
the boundary values away from <9ct+ U der_ exist and we have 



(5.21) h T (X n ,n) — h T (X l ,n), Aecr + Ua_. 
By property I, (b), 

(5.22) h =F = T± 1 (i2 ± v±+^) = ^ + ^eR, Xeo-^, 



from which 

(5.23) h T (X n , n) — h T (X l , n), Agct^, 

follows. Combining (|OTj) and (jST^j) yields ((5TT4")) . 

(hi). Since the last two terms in (|5.12p are 0(z~ 1 ), the asymptotic behavior 
follows from (|3.5p and II, (c). 

(iv). From (|5TT31) . ([3~5|) . and (|3~Tg|) we obtain 



„ ±(M ^ ±W) . EMpi) . _J_ . ± „ Wi Ae , ± . 

Hence equality holds for all z G C by analytic continuation. □ 
Corollary 5.6. The functions hzf(z,n) admit analytic extensions to C\ a^. 

Proof. Property (i) of Lemma [5Tl holds for z E C \ a. Relation (|5. 14|) implies that 
hzf: have no jumps across z € wk(cr± ). To finish the proof we need to show that 

the possible remaining singularities at E G dcr^ n da are removable. This follows 
from (cf. (15161) 1 

(5.24) Ct M = T =S=4 (z, n) 

which shows C^(z,n) — 0((z — E)^ 1 / 2 ) and hence h T (z,n) = 0((z — i?)" 1 / 2 ) for 
E G <r£ X) n 9cr. 

However, let us emphasize at this point that the behavior of h±{z,n) at the 
remaining edges is a more subtle question to be discussed later. □ 

Eliminating ip± from 



J R ± (A)V±(A,n) + V±(A,n) = h T (X,n)T±(X) 



fl±(A)V±(A,n)+V±(A,n) = /i T (A,n)T±(A) 
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yields 

V±(A, n) (1 - |i?±(A)| 2 ) = h T (X,n)T ± (X) -R^X)h T (X, n) T± (A). 
We apply I, (c), II, and the consistency condition I, (d) to obtain 



(2) 



T T (A)^ ± (A,n) = h T (X,n) R±i ^ ± W h T (X,n) 

(5.25) = h T (\,n) + R T (\)h T (\,n), X G a 

This equation together with (|5 . 13f> gives us a system from which we can eliminate 
the reflection coefficients R±. We obtain 

(5.26) 

r±(A)(V±(A)V T (A)-h±(A)fc T (A)) =V±(A)/i±(A)-<MA)MA), Aei^K 

Now introduce the function 

(5.27) G(z) := G(z, n) = ^(^")^-(^^-M^^-(^") 

which is well defined in the domain z G C \ (a U a d U M+ U M_). By (fST2"6")) and 

(5.28) G(A) = (^±(A)MAj-^(A)^±W)p±(A) J A e 

so we need to study the properties of G(z, n) as a function of z. Our aim is to prove 
that G(z,n) = 0, which will follow from the next lemma. 

Lemma 5.7. The function G(z,n), defined by f|5 . 2T[) . has the following properties. 

(i) G(A U , n) = G(A\ n) G K /or A G R \ (9cr_ U cV+ U cr rf ). 

(ii) -ft /ias removable singularities at the points da- U da+ U <7d, where ad '■= 
{Ai, A p }. 

Proof, (i). We can rewrite G(z,n) as 

ip-\-(z, n)ip-(z, n) — h+(z, n)h-(z, n) 



(5.29) G(z,n) 



W(z) 



where h±(z, n) — 5±(z)h±{z, n) as usual. The numerator is bounded near the points 
under consideration and the denominator does not vanish there. Thus G(z,n) has 
no singularities at the points (M+ U M_) \ ad- 

Furthermore, by Lemma 15.41 II, (a), and Lemma 13. II we know that G(z,n) has 
continuous limiting values on the sets cr_ and <r+, except possibly at the edges, and 
satisfies 



G(A u ,n) = G(A 1 ,n), Aeo- + Uo-_. 

Hence, if we can show that these limits are real, they will be equal and G(z, n) will 
extend to a meromorphic function on C, that is, (i) holds. To this aim we first 
observe that (|5.14p . (|5.28p . and Lemma l3~Tl imply 

(5.30) G(A U , n) = G{X\ n) G R, A G int(a (2) ). 
Thus, it remains to prove 

(5.31) G(A u ,n) = G{X\n) G K for A G mt(cr^) U in.t(<7+^). 
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Let us show that G(A, n) has no jump on the set int(<yi 1 ' > )Uint(<7+ 1 ' ) ). We abbreviate 



(5.32) 



[G] := G(A) - G(A) = 







h + h- 


w 




w 



A £ 



and drop some dependencies until the end of this lemma for notational simplicity. 

Let A £ int(<4 1),u ), then ip±,h± £ R and T T = -(Wp^)' 1 . By 1(05)1 . (I), (b), 
and (|5.13[) we obtain for A £ int(<rl 1 '') 



(5.33) 



■0± 



(l),u. 



Since p± £ R for A e int(a^ hn ), ([OJD implies 



The only non-real summand in (|5.12[) is the Cauchy-type integral. The Plemelj 
formula applied to this integral gives 

[h T T ± ]=-p T ^ ± \T T \ 2 W(^,^ ' " ' 



P± 



and by (|5.33p we get 
(5.34) 



h + h- 




~1p+1p-~ 


w 




w 



p T ^ ± h ± \T T \ 2 , AGint(4 X) ) 



bince W ^ for A £ in^ai^), the function 



p T ip±h ± \T T 



2 _ 5^ ^ ± l ± 



Pt \W\ 2 

is bounded on the set under consideration. Finally, (|5.34[) and (|5.32[) imply (|5.31[) . 

(ii). Now we prove that the function G{z,n) has removable singularities at the 
points <9cr_ U da + U ad- We divide this set into four subsets 

(5.35) (if = <9ct (2) n int((T T ), n 2 = <9ct (2) n da, n± 



flf = da v ±>nda ± , 



to 4 



Since all singularities of G are at most isolated poles, it is sufficient to show that 

(5.36) G{z) =o({z-E)~ 1 ) 

from some direction in the complex plane. 

to\: Consider E £ to* (the case E £ f2j~ being completely analogous). We will 
study limA^E G(A, n) as A € int((7^ 2 ^) using (|5.28[) with the "— " sign. Note that 
ip- = 0(1), p- = 0(1), and W(E) ^ 0. Moreover, we obtain from Lemma [3.11 
respectively II that 



which shows 



0(1), 



o 



fX^E 



E i M+, 
E £ M+, 







T+(X) 



X-E 

0(1), 



E 4 M. 



MA) 



^ + (A) + i? + (A)^+(A) 



r+(A) 

for A £ (j( 2 ) . Inserting this into (|5.28|) shows G(A, n 
case E £to\. 



O 



1 



\/A" 
0(- 



E£M + , 



and finishes the 
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VL 2 : For E £ da {2) n da, we use ([S~28| and take the limit A -> J5 from cr( 2 ). First 
of all, observe that 

The case E G er„ is evident. If ^ cr„ then (f^TTT]) and (pT23|l yield 

Therefore, both for virtual and non-virtual levels the estimate 

(5.37) 6- (R-ip- +i>J) W = o(l), Eeda-, 

is valid. Inserting (|5.13[) into the summand ip + h + p + of (]5 . 28[) (for the second 
summand we use an analogous approach) we obtain (recall l|2.2[) ) 

(5.38) =|+^L(i2_V-+?-)W- 
Combining the estimate 

^+ f 1 



with (|5.37p we have G(z) — o((z — E)^ 1 ) as desired. 

^3: Suppose that E S 9e/_ 1 ' ) n9<7-_ (the case E € da^ nda + is again analogous). 
Now we cannot use (15.281) . so we proceed directly from formula (|5.27p estimating 
'+'/'- anr i h + h - 

the set int (c_ ). By Lemma [37X1 and (|3 . 22[) we have 



the summands ^tjf and h+ J [ t separately. We investigate the limit as A — > E from 



(5.39) ^_ = ^_ = Q ( I 



W w VvA - 

hence the first summand has the desired behavior. To estimate the second sum- 
mand, we split the function h—(X,n) according to 

h-(\,n) = hi(X,n) + /i 2 (A,n), 

where 

h 1 (X,n)=W+ n _ 1 (C-(X,-),d^(X,n,-)), h 2 (X,n) = h-(X,n) - Ai(A,n), 

(5.40) d_(A,n,.) :=y u — <&>-(0- 

It follows from the proof of Lemma T5. 41 that /12(A) = 0(£_(A)) for A — ► -B. Recall 
that at the point under consideration singularities -E € {m^", • ■ ■ , A 1 ^} U might 
occur (in the case da^Dda one can have E 6 M + \JM + and in the case da^ C\da^ 
one can have E e -M+). Introduce 

(5.41) #(*,n) :=$ + (z)$+(a! s n) 
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and recall that (12.15P implies 



(5.42) <f>+(z 7 n)-0+(E,n) = O(V^E). 
Then (see ([2J2]) and (f27T3|) > we have 

(5.43) b±^ =0 (h^l) =0 (^±^l) =o( h+ ^ 



W V WT+ I V P+ / P. 



Now we distinguish two cases: (a) E G da^ n 9o+ and (b) E S <9eri 1 ' ) n da. 
Case (a). By (|5.13l) and (|5.37|) we have 

(5 44) 5 + / 1+ = ^- V '- + ^-^ + ^-(«-^-+^-) Y 1 



T- P- yy/z^Ej' 

therefore 

(5.45) h+{ 5-} X) =o( 1 ^ 



w(A) VVa~^b;p+(a) 

As a consequence of |r- = 0(^j===) we obtain 

(5.46, ^-(rb). 

Next, we have to estimate 

By (|5.42p we can represent (|5.45|) as 

( 5 48) ^(A)C-(A) _ / 1 \ / ^+(g) 

(5 ' 48) w(A) - 1 7x=^ + { 

Then (|5T47)) implies 

h + (X,n)hi(X,n) f 1 



, x W{X) VVA^B 

(5.49) 



0(d_(A,n)) + 0(d_(A,n-l)) 



Q«n-1 (<#(£), rf-(A))) \ 

To estimate d_ in the first two summands we distinguish between the resonance 
case, E G <7„, and non-resonance, i? ^ a v . First let .E ^ a v , that is, VF(E') 7^ 0. 
From (|5.19p and (|5.20p we see that the integrand is bounded as A — > E ^ <7„, then 
d_(A) = 0(l)by[T3. 

If E e cr t ,, then (|*H2]) (see also (EH])) yields 



|r_(0lV(0<M£, •)<(£, = o 

and [171 Eq. (29.8)] implies 



1 
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For the estimate of the last summand in (|5.49|) we use (|5.19|) and (|5.40[) to represent 
the integrand in W^ n _ 1 ($+(E),d-(Afj as 

\T-(t)\*P-(M+&n)Kn-i 

r- ^ 



It follows from (f27T5|) and (|OT|) that 



(j&2^W<o.tf(*>)). 



Wi,-iW(f),«E)) = O(V^TE), 

which implies together with (I3.22|) the boundedness of the integrand near E. Thus, 

(5.51) W+ n _ 1 (4,+ (E),d.(X)) =0(1), 

and combining ([5T6|) . (|5~49]) . ([5~5Tjl . and ([5~50| finishes case (a). 

Case (b). Now we do not have estimate (|5.37[) (cf. Ill, (b)) at our disposal, but 
we can proceed as in (|5.43p , (I5.44[) since P+(E) ^ and arrive at 

(5JB, ^ =0 (M + ) = 0f«^±±t)U f * 



This estimate is sufficient to conclude that ()5.46j) is valid in case (b) as well. For 
hi, we use the following estimate (cf. (|5.50|) and (15.52| ) instead of (I5.47|) : 

h + hi (h£\ 0{d _ ) = ( W \J 1 



w v w j v ' \^\~e) kvx^e 

Combining this with (|5.46|) finishes case (b). 

SI4: Finally we have to show that the singularities of G(z, n) at the points of the 
discrete spectrum are removable. Since W(z) has simple zeros at z = Afc, it suffices 
by (j5~2"5)) to show that 

(5.53) h + (Xk,n)h-(Xk,n) = ^-(\ kl n)^ + (\ k ,n). 

By Lemma l5.4[ the functions hzp = S^h^ given in (|5.12p are continuous at the points 
M±. Since (S T T± x )(\ k ) = and (S^T^ 1 ^f)(X k ) = 0, only the last summand in 
(I5.12j) is non-zero. We compute the limit of this summand as A — ► X k using (|3.19D . 

(5.54) MA fe ) = -7 ± ,fe^ ± (A fe )^^, 

and apply (|3~20|) to obtain (|5~53)) . □ 
The identity G(z,n) = implies 

(5.55) Vm-( z j n)"0- {z, n) — h+(z, n)h-(z, n) = 0, Vn £ Z. 
For z^cowc obtain by (|2.16|) and (|5.9p 



1> + (z, n)1>-{z, n) = K+(n, n)K_ (n, n) J] * -^(1 + o{\)), 

j=o a <? U) 
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Formulas 15. 15[) and ()3.23|) imply 

T+(oo) 2 K + (n,n)K-(n,n) jJ- Q a^{j) 



and by (|5.55[> . 

fJo a i 0) T +(°°) 
The value on the left hand side does not depend on n, so using (|5.7[) we conclude 

(5.56) a +( n ) = a-{n) = a(n), Vn € Z. 

It remains to prove b+(ri) — b—(n). If we eliminate the reflection coefficient R± 
from (|5.13p at n and (|5.25|) at n + 1 we obtain 

^+(A,n)^_(A,n + l) - /i+(A,n + l)/i_(A,n) 

(5.57) W(A) 

= p+ (A) (fc± (A, n + l)ip± (A, n) - V± (A, n)ft± (A, n + 1)) , A e a* 2 ^' 1 . 

Proceeding as for G(A, n) in Lemma [5.7l we can show that that the function G\ (z, n) 
is holomorphic in C. From (pHS) . (|5J| . (pT2"3"]) . (|2~T6| . ([5756]) . and the Liouville 
theorem we conclude that 

ip+(z, n)ip-(z, n + 1) - h + (z, n + l)/i_(z, n) . 

W{zr~ = ~ 1/a{n) - 

We compute the asymptotics of 

W(z, n) :— a(n) (i()+(z, n)ip-(z, n + 1) — h + (z, n + l)h-(z, n)) = — W(z) 
as z — * oo and obtain (compare (|3.5[) 1 

(5.58) = W(z,n) - W(z,n-1) = (b+(n) - b-(n))K + (0,0)K-(0,0). 

This implies in particular b+(n) — b-(n) = 6(n), hence the proof of Theorem l5.3l is 
finished. 
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